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Abstract 

Prudent walks are self-avoiding walks which cannot step towards an already occupied 
vertex. We introduce a new model of adsorbing prudent walks on the square lattice, which 
start on an impenetrable surface and accrue a fugacity a with each step along the surface. 
These are different to other exactly solved models of polymer adsorption, like Dyck paths, 
Motzkin paths and partially-directed walks, in that they are not trivially directed - they are 
able to step in all lattice directions. We calculate the generating functions, free energies and 
surface densities for this model and observe a first-order adsorption transition at the critical 
value of the surface interaction. 


1 Introduction 

Lattice walk models have long been considered as a canonical model for long-chain polymers 
in solution [21]. One popular lattice walk model in the polymer community is that of self¬ 
avoiding walks (SAWs). The standard model associates a weight (or fugacity) z with each step, 
or monomer. If Cn is the number of n-step walks (equivalent up to translation), the length 
generating function is 

F{z) = J2cnZ'^- 

n>0 

It has long been known that the limit 

log y = lim — log Cn 

n^oo Ti 

exists, where the growth constant y depends on the lattice being considered. For most lattices in 

2 or higher dimensions, only numerical estimates for y are known [15,16]. A notable exception 
is the 2-dimensional honeycomb lattice, for which it has recently been proven by Smirnov and 
Duminil-Copin [8] that y = \/2 + \t2. 

To model the behaviour of a polymer interacting with an impenetrable surface, we consider 
walks restricted to a half-space. Although this restriction limits the number of walk configura¬ 
tions, it has no effect on the exponential growth of the number of walks and hence the growth 
constant [27]. To model the monomer-surface interactions of an adsorbing polymer, we associate 
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a fugacity e“ with every edge of the walk which lies in the boundary of the half-spaced If we 
define as the number of n-step half-space walks beginning on the surface with v edges 

along the surface, then the partition function is given by 

n 

Zn{a) = ^c+(zy)e^“, 

u=0 

where a = —e/ksT. Here e is the energy associated with a surface visit, T is the absolute 
temperature and ks is Boltzmann’s constant. For sufficiently large values of a, configurations 
with many visits to the surface dominate the partition function and the walk is said to be in 
an adsorbed state; otherwise, the loss in configurational entropy dominates, and the walk is 
repelled by the surface and said to be in a desorbed state. 

The free energy of the system is 

/(a) = lim - log Zn{a), 

n^oo Ti 

which has been shown [11] to exist for all a < oo. When a < 0 the free energy is independent 
of a and given by log/x [27]. For a > 0, it has been shown that in two dimensions 

/(a) > max{log/x, a}. 


This implies the existence of a critical value etc with 0 < etc < log/x where the free energy /(et) 
is non-analytic. We identify this point of non-analyticity as the location of the adsorption phase 
transition: for a < Oc, the polymer is desorbed, while for et > etc the polymer is adsorbed. 

Closely related to the free energy is the mean density of surface visits; in the limit of infinitely 
long polymers this is given by 


p{a) = 


dfja) 

da 


( 1 ) 


This variable acts as an order parameter for the system, and signals the onset of a phase 
transition. If a < Oc then p{a) = 0, while p{a) > 0 for a > Uc- If the adsorption transition 
is second-order then p is continuous for all a, while a first-order transition is manifested by a 
jump discontinuity in p{a) at Uc- 

Instead of working directly with the partition functions, Zn, we will consider the generating 
function Z{z,a) = Yln^n{oi)z'^. A central tenet of analytic combinatorics [9] is that the dom¬ 
inant singularities (the points of non-analyticity closest to the origin) of a generating function 
determine the asymptotic behaviour of its coefficients. Specifically, if for a given a the dominant 
singularity of Z{z, a) occurs at z = Zc(a), then 


f{a)= lim -logZ„(a) =-logZc(a;). (2) 

n—>-oo Ti 

For many models (including the one considered in this paper) it is much easier to determine 
the generating function than to compute the individual partition functions Z^ for each n. But 
since (2) enables us to determine the free energy /(a) of a model directly from its generating 
function, we need never consider the Zn anyway. 

For ease of notation we will frequently use the surface fugacity a = e“, and consider the free 
energy, etc. as functions of a. The above arguments are still valid, with the critical fugacity 
occuring at Oc = e“‘'. 

^One can instead associate the fugacity with vertices in the surface, to obtain a slightly different model. 
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Figure 1: Examples of directed polymer adsorption models, with vertex weights indicated: (a) 
a Dyck path, (b) a Motzkin path and (c) a partially-directed walk. 

1.1 Exactly solved models 

At the present time, a closed form expression for the generating or partition function of SAWs 
is still unattainable. Nnmerical experiments, using series analysis [5, 10] or Monte Carlo simu¬ 
lations [4, 14], can be used to estimate non-rigorous values of // and Oc, for example. 

Exactly solvable models can be obtained by placing a restriction on the walk models used to 
study the polymer system. Possibly the simplest restriction to nse is directedness - forbidding 
steps in certain lattice directions. Below, we summarize the three most commonly-used directed 
models. 

• Dyck/ballot paths. [12] These are generated on the (45° rotated) square lattice by 
forbidding NW and SW steps. Since there are no steps parallel to the surface, only 
vertex-weighted models are possible. In the desorbed phase, these walks have growth 
constant /i = 2. The adsorption phase transition occurs at Oc = 2. 

• Motzkin paths. [13] These are generated on the triangular lattice by forbidding NW, W 
and SW steps. The growth constant is fj, = 3, and the phase transition occurs at Oc = 2 
for edge weights and Oc = 3/2 for vertex weights. 

• Partially-directed walks. [25] These are generated on the square lattice by forbidding 
W steps. The growth constant here is ^ = 1 -|- \/2, and the phase transition occurs at 
Oc = (2 -|- \pi)l2 for edge weights and Oc = (1 -|- \/2)(\/5 — l)/2 for vertex weights. 

All of these directed models undergo second-order adsorption transitions at their critical 
valnes of the surface fugacity. These models are relatively easy to solve, and can easily be 
adapted to model other polymeric objects and phenomena, such as copolymers and polymer 
collapse [3,18,23]. The directedness restrictions are, however, very strong, and result in walks 
which have considerably less freedom than general SAWs. 

Here we introdnce a model of polymer adsorption on the square lattices which does not 
require a directedness restriction - the walks are able to step in all lattice directions. Instead of 
being directed, the walks must be prudent - they can never take a step towards a previously- 
visited lattice point. The model of partially-directed walks can be viewed as a subset of this 
prudent walk model. 

Using the prudence restriction we recursively generate the walks by nsing a set of generating 
functions that keep track of certain geometric measnrements. These recurrences can be written 
in the form of fnnctional equations, which can be solved via the kernel method and some of its 
generalisations. 

We analyze the generating functions to determine the free energy /(a) and the surface 
density p{a) for prudent walks. These calculations are carried out for walks where either one 
or both endpoints of the walk he in the surface and we refer to these two snbsets as tails and 
loops, respectively. We hnd that both prudent tails and loops undergo a first-order adsorption 
transition at a(, = 2 and ~ 1.82476, respectively. 
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Figure 2: Prudent walks on the square lattice; (a) 2-sided, (b) 3-sided and (c) 4-sided (unre¬ 
stricted) . 


1.2 Prudent walks 

Prudent walks on the square lattice were introduced by Prea [24] and later studied by Duchi 
[7] and Bousquet-Melou [2]. These walks have the useful property that their endpoint always 
lies on the boundary of their bounding box - the smallest lattice rectangle which contains the 
walk. Using this property, Duchi dehned four classes of increasing generality: 

• 1-sided prudent walks must always end on the E side of their bounding box, 

• 2-sided prudent walks must always end on the N or E sides, 

• 3-sided prudent walks must always end on the N, E or W sides (and may not step between 
the SE and SW corners of the box when the box has width one), and 

• 4-sided (or general) prudent walks may end anywhere on their box. 

The class of 1-sided prudent walks are equivalent to partially-directed walks, a model of 
polymer adsorption that has been studied for some time (for an early example see [25]). Duchi 
and Bousquet-Melou solved the generating functions of 2- and 3-sided prudent walks (without 
a boundary) respectively. The generating function of 2-sided walks is algebraic, while that 
of 3-sided walks is non-D-finite. Bousquet-Melou also wrote down a functional equation for 
the generating function of general (4-sided) prudent walks, but no solution for the generating 
function has yet been found. 

2 Generating functions 

We are able to solve the functional equation for 2-sided prudent walks above an impenetrable 
surface, however the addition of the surface and the corresponding interaction fugacity a makes 
the generating function significantly more complex than that of the original (surface-free) model. 
Determining the critical critical behaviour of the adsorption model requires a great deal more 
effort. 

We have already noted that 1-sided prudent walks are equivalent to partially-directed walks. 
Eor the 3-sided model, we are able to write down functional equations, but the number of 
catalytic variables is too great to obtain a closed form solution. There are no 4-sided walks which 
are not also 3-sided walks. To see this, note that the S side of the bounding box necessarily lies 
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Figure 3: 2-sided prudent walks above an impenetrable surface, (a) ending on the right of the 
box and (b) ending on the top, with the distances i and j indicated. 


on the impenetrable surface. So any step taking the walk from the E or W sides onto the S side 
must be in the direction of the origin, and thus cannot possibly be prudent. 

For the 2-sided purdent walk model, we begin by defining the generating functions 

R{u,v) := R{z-,u,v,a) = Rn,ij,yZ^u''v^a'' 

T{u,v) := T{z;u,v;a) = Tn^ij^yZ^u^v^ 


where Rn^ij^u (resp. is the number of n-step 2-sided prudent walks which start on a 

horizontal impenetrable surface and end on the right (resp. top) side of their bounding box, 
with a distance i from the endpoint to the top (resp. right) of the box, a distance j from 
the endpoint to the surface, and u steps along the surface. We write R{u,v) and T{u,v) for 
shorthand as u and v are the catalytic variables - the ones used in the iterative construction of 
the walks. The introduction of the surface breaks the symmetry between walks ending on the 
top and right of the box, forcing us to use two generating functions. 


Lemma 1. The generating functions R{u,v) and T(u,v) satisfy the functional relations 

1 


L{u, v)T{u, v) = 


z^v 


1 — zua u — z 


T{z, v) -|- zR{zv, v) — z{l — a)R{zv, 0) 


z'^v 


ziii 1 — d) 

M{u, v)R{u, u) = 1 -|- zvT{z, v) -— R(zv, v) - R(u, zu) - -R{u, 0) 


z‘^u 


u — zv 


V — zu 


+ 


u — zv 
z‘^v{l — a) 


u — zv 


R{zv, 0) 


where 


L{u, v) 
M{u, v) 


L{z; u,v) = 1 — 
M{z] u^v) = 1 — 


zuv{l — z'^) 

{u — z){l — zu) 
zuv{l — z^) 

{v — zu){u — zv) 


(3) 

(4) 


Proof. We iteratively construct walks by considering their last inflating step - the last step 
which moved either the top or right side of the bounding box. We begin with the construction 
for walks ending at the top of their box. There are three possibilities: 


• Case 1. Neither the top nor right sides of the box have ever moved. Then the walk is 
empty or consists only of W steps. The generating function for such walks is 

1 I I 2 2 2 I 1 

1 -|- zua + z u a -|- ... =-. 

1 — zua 
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• Case 2. The right side was the last to move. Then immediately prior to this, the walk 
ended on the right side of its box, and it then took an E step. After this it must have taken 
N steps to reach the top of the box, but no further. There are, however, two sub-cases 
to consider, depending on whether that E step was along the surface. If it was then the 
walk gains an a weight; otherwise it does not. 

Case 2a. The inflating E step was along the surface. The generating function for 
such walks is 

= zaR{zv,0). 


Case 2b. The inflating E step was not along the surface. The generating function 
for such walks is then 

a’^ = zR{zv,v) — zR{zv,0). 

nA.jM 

j>0 


• Case 3. The top side was the last to move. Then immediately prior to this, the walk 
ended on the top side of its box, and it then took a N step. After this, it may have taken 
E steps as far as the right side of the box, or any number of W steps. The generating 
function for such walks is then 

( i oo \ 

fc =0 1=1 / 


(where the sum over k counts the walks which go E after the N step, and the sum over I 
counts those which go W after the N step), 


E 




zv 


- {uT{u, v) - zT{z, v)) + / 

u — z 1 — zu 


- 1 “- 

u — z 1 — zu 

T{u,v). 


Every walk counted by T{u,v) falls into exactly one of Cases 1-3. We thus have 


T{u,v) = 


1 


1 — zua 


+ zaR{zv, 0) -|- zR{zv, v) — zR{zv, 0) 




zv 


u — z 


{uT{u,v) — zT{z,v)) 


z^uv 
1 — zu 


T{u,v). 


Rearranging this gives (3). 

A similar consideration can be made for walks ending on the right side of their box. The 
same three cases apply. 

• Case 1. Neither the top nor right sides of the box have moved. The only such walk is 
the empty walk, with generating function 1. 

• Case 2. The right side of the box was the last to move. Then the last inflating step was 
E. After this, the walk may take S steps as far as the surface or N steps as far as the top 
of the box. 
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Case 2a. The last inflating step was along the surface, and the walk gained an a 
weight. The generating function for such walks is 

ZOj 

- iuRiu, 0) — zvR{zv, 0)). 

u — zv 

Case 2b. The last inflating step was not along the surface. The generating function 
for such walks is 

z z 

{uR{u, v) — zvR{zv, u))- {uR{u, 0) — zvR{zv, 0)) 


u — zv 


u — zv 


z^u 

H- {R{u, v) — R{u, zu)). 

V — zu 

• Case 3. The top side of the box was the last to move. Then the last inflating step was N. 
After this, the walk stepped E as far as the right side of the box. The generating function 
for such walks is 

zvT{z, v). 

Summing all these contributions, we obtain 

ZOj Z 

R{u, u) = 1 H - {uR{u, 0) — zvR{zv, 0)) H - {uR{u, v) — zvR{zv, v)) 


u — zv 


u — zv 


{uR{u, 0) — zvR{zv, 0)) + 


z‘^u 


{R{u, v) — R{u, zu)) + zvT{z, v) 


u — zv V — zu 

which is equivalent to (4). □ 

Lemma 2. The generating functions T(u,v) and R{u,v) satisfy the functional equation 

(5) 


where 


M{u, v)R{u, v) = A{u, v) + B{u, v)R{u, 0) + C{u, v)T{z, v) 
1 + z"^ — zv + z^v — 


A(u) := A(z; v) = 

A{u, v) := A{z‘, u, V, a) = 


B{u, v) := B{z; u, v;a) = — 

C{u, v) := C{z] u,v) = — 
Proof. Setting u = 0 in (4) yields 


2z 

v{u — z^u — zuX{v)a + z^vX{v)a) 

{u — zv)(v — zu)(l — zX(v)a) 
zu(u + V — zv — z'^v — va + z‘^va) 
{u — zv)(v — zu) 
zv(zu — uX(v) + zvX(v)) 


(X(v) — z)(u — zv) 


{\ + z — za)R{u, 0) = 1 + zR{u, zu), 

which can be used to eliminate R{u,zu) from (4). Meanwhile, X{v) satisfies the equation 
L{X{v),v) = 0 (as does another function of v, but A is the only root which is a power series 
in z and thus the only one which will eventually give a well-defined solution), so substituting 
u = X{v) into (3) cancels the LHS. The resulting equation can be written as 

z‘^v{l — a) 


z‘^v 


zv 




u — zv 


u — zv 


{u — zv){l — zX{v)a) {u — zv){X{v) — z) 


■R{zv, v) H- -R{zv, 0) = ^ ^ ^ ^ ^ ^ - ^T{z, v). 

( 6 ) 

This can be used to eliminate R{zv,v) and R{zv,Q) from (4). Simple manipulation gives (5) in 
the stated form. □ 
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( 7 ) 


Proposition 1. The form of R{u,0) is given by 


oo n—1 

R{u,0) = 

n=0 k=0 

when eonsidered as a formal power series in z, u, a, where 


A:=A(1) = 


1 - Z + + ^3 _ Vl - 2z - z2 _ ^4 2z^ + 

1 (1 - z2)(l - A2)(l - zail + z)X{qA) + zaX{qAf) 


niq) := n{z]q;a) = 

Z{q) := Z{z- q- a) = _ ^a)( 1 - (1 - z - z^ + A)A(gA)) 

i7(a) := J'(z; a) = 1 + A — zA — z^A — Aa + z^Aa 


J{a) (1 — zA)(l — zaA(gA))(l — [1 — z — z"^ + A)X{qA)) 

1 (A — z)(l — z — — a + z^a + A)(z — (1 — zA)A(gA)) 


Proof. Equation (5) is susceptible to the iterative kernel method. The kernel M is cancelled at 
{u,v) = {u,uA) (and since M{u,v) = M{v,u), it also disappears at {u,v) = {vA,v)), and so by 
substituting the pairs 

{u,v) = {u,uA) and {u,v) = {uA'^ ,uA), 

we obtain two equations in R{u,0), R{uA'^,0) and T{z,uA). Eliminating T{z,uA) gives 

R{u,0) = T-L^u) + Z{u)R{uA‘^ ,0) ( 8 ) 


where 

, , A{u,uA) C{u,uA)A{uA‘^,uA) 

B{u,uA) B{u,uA)C{uA'^,uA) 

. , C{u,uA)B{uA^,uA) 

B{u,uA)C{uA‘^ ,uA) 

After simplification R and Z take the form given in the Proposition. 

Iterating (8) gives (7), provided that everything converges as a formal power series. But now 
R{u) is a power series in z with coefficients in Z[ri, o] of the form 1 + za + O(z^), and likewise 
Z{u) is a power series of the form z^(l — a — u + ua) + 0(z®). As such, the sum converges as a 
formal power series in z with coefficients in Z[ri, a]. □ 


Theorem 1. The generating functions R{u,v) and T{u,v) have the solutions 
R{u, v) = [A{u, v) + B{u, v)R{u, 0) + C{u, v)T{z, u)] 

1 

u — z 

where 

T{z, v) = [A{vA, v) + B{vA, v)R{vA, 0)] 

and R{u, 0), R{vA, 0) are given by (7). The overall generating function for 2-sided prudent walks 
above an impenetrable surface is given by 

W{u,v) := W{z-,u,v,a) = Wn,i,j,vZ^u''v^a'^ 

= i?(u, v) + T{u, v) — r(0, v) 

where Wn,i,j,v is the number ofn-step 2-sided prudent walks which end a distance i from the NE 
comer of their box and a distance j above the surface, with v steps along the surface. 



T{u,v) = 


L{u, v) 


1 


1 


1 — zua 1 — zA(u)a 


+ z^v 


1 


X{v) — z 
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Figure 4: 3-sided prudent walks above an impenetrable surface, (a) ending on the right of the 
box and (b) ending on the top, with the distances i, j and k indicated. 


Proof. Substituting {u, v) = {vA, v) into (5) cancels the LHS, and rearranging gives the stated 
form of T{z, v). Since vA is a power series in z with coefficients in Z[u], R(vA, 0) is well-defined 
as a formal power series. Rearranging (5) then gives R{u,v). 

Equation (6) can be rewritten as 

— 1 z^v 

zR{zv, v) - z(l - a)R(zv, 0) = -- -h —rw- T(z, v), 

1 — zX{v)a X[v) — z 

and this can be used to eliminate R{zv,v) and R{zv,0) from (3). Rearranging the resulting 
equation gives the stated expression for T{u, v). 

The overall generating function W{u,v) is then found by adding the generating functions 
for walks ending on the right and the top of the box, and we subtract r(0, v) (or equivalently, 
R{0,v)) because walks ending at the NE corner have been counted twice. □ 


We present here the functional equations for general (3- or 4-sided) prudent walks on the 
square lattice, though we are unable to solve them. Define the generating functions 


R*{u,V, w) 
T*{u,v,w) 


R*{z] u, V, w] a) 
T*{z] u, V, w] a) 


n,i,j,k,u 


where R* ijku counts n-step walks ending on the right of their box and counts those 

ending at the top of their box. In both cases i is the distance from the endpoint to the NE 
corner of the box, j is the distance from the endpoint to the surface, k is the distance from the 
endpoint to the W side of the box, and v is the number of occupied edges along the surface. 

Lemma 3. The generating functions R*{u,v,w) and T*{u,v,w) satisfy the equations 


L*{u, V, w)T*{u, u, tc) = 1 -|- zwR*{zv, v, w) + zw{a — l)R*{zv, 0, w) -|- zuR*{zv, v, u) 

7‘^nui z'^un (H) 

+zu{a — l)R*{zv, 0, u) - T*{zw, v, w) - T*{u, v, zu) 

u — zw w — zu 


z^vw 


L*{u, w, v)R*{u, v,w) = 1 - R*{zv, v, w) — 


z^uw 


u — zv 


V — zu 

v2 


N zuwia — 1) , 

R {u, zu, w) -|- -R {u, 0, w) 


u — zv 


z^vw{a — 1) 


u — zv 


R*{zv, 0, w) -|- zvT*{zw, V, w) 

( 10 ) 


where 


L*{u, V, w) := L*{z] u,v,w) = l — 


zuvw{l — z'^) 

{u — zw){w — zu) 
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We omit the proof as the details follow in much the same way as Lemma 1. Here, inflating 
steps can be N, W or E; inflating steps to the N are attached to walks counted by T*, and 
inflating steps to the W and E are attached to walks counted by R* (symmetry means R* also 
counts walks ending on the left side of the bounding box). 

Though the equations are somewhat more complicated, the problem we face here is similar 
to the difficulty of solving general prudent walks without a boundary [2] - the presence of 
three catalytic variables seems to render the kernel method unable to produce a closed form 
expression. 


3 Singularities and asymptotic behaviour 

While we are quite certain of the location of the dominant singularity of W{z\l,l]a) for all 
a > 0, the complicated nature of the generating functions has prevented us from obtaining a 
rigorous proof. The results of this section are thus stated as conjectures. They are followed at 
the end of the section by an outline of our reasoning. 

Conjecture 1. For a given a > 0, the dominant singularity ofW{z; 1, l;a) is located at 

tf \ ^ 2 

z (a) = { ! 

a>2, 

where z\ k, 0.403032 is a root of 1 — 2z — 2z'^ + 2z^ = 0, and z^ia) is a root of 


1 — a — a(l — a)z + az^ + a(l — a)z^ = 0. (11) 

We thus observe a first-order adsorption phase transition occuring at a = 2, with a crossover 
exponent of cj) = 1. 

See Figure 5 for plots of the dominant singularity and surface density, as calculated using (1) 
and (2). It is of interest to note that the equation which determines the free energy in the 
adsorbed state, namely (11), is the same for prudent tails and partially directed walks [28]. 

Series analysis of W{z;l, 1; a) for a variety of a values confirms the validity of this conjecture 
- see Figure 6 for a comparison of the free energy and corresponding numerical estimates. We 
point out that the dominant singularity in the desorbed phase, z\ ~ 0.403032, is the same as 
for 2-sided prudent walks without a boundary [2] . 

Conjecture 2. For a given a > 0, the dominant singularity ofW{z; 1,0; a) is located at 

f, ^ {z{ a < PS 1.82476 

[z^[a) a > al, 

where z{ ^ 0.412095 is a root of 

l-3z-z‘^+ 6z^ - 7z^ -z^ F 3z® + = 0, 

is a root of 

1 - 7a + 450^ - 143a^ + 2770^^ - 346a® -h 285a® - 155a'^ + 54a® - lla^ + a^° = 0, 

and Z 2 {a) = z\{a) as defined in Conjecture 1. This implies that in the absorbed state, the free 
energy for prudent tails, prudent loops and partially directed walks is determined by the same 
equation, namely (11). We thus observe a first-order adsorption transition occurring at a = 
and a crossover exponent of (j) = 1. 
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Figure 5: Top: The (conjectured) free energy of 2-sided prudent walks (solid) and loops (dashed) 
as functions of the surface fugacity a = logo. Bottom: The (conjectured) density of edges in 
the surface for 2-sided prudent walks (solid) and loops (dashed). Note the discontinuity in the 
surface density for both models, corresponding to a first-order adsorption transition. 



/ 


(See Figure 5 for plots of the free energy and surface density.) 

The dominant singularity of W{z; 1,1; a) in the desorbed phase, namely z\, appears as a pole 
in the T{z; z, 1; a) term, which appears in both R{z] 1,1; a) and T(z; 1,1; a). More specifically, 
z{ is a root of C(z; A, 1) = 0. The dominant singularity in the adsorbed phase appears as a 
singularity of R{z; 1,0; a) and R{z; A, 0; a). In fact, z^ia) is a root (in the variable z) of J'{a) = 0 
(as defined in Proposition 1), which in turn makes it a root of B{z; 1, A; a) = 0 (and in fact a 
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Figure 6: A plot of the conjectured free energy f{a) (the red line) along with numerical estimates 
(the blue dots) computed using series analysis of walks up to length 400. 


root of i?(z; a) = 0 for n > 0), and thus a pole of 7i{z; a) and Z{z\A^]a) for all 

n > 0. 

The trouble here lies not in locating the singularity given in Conjecture 1, but rather in 
proving that kF(z;l,l;a) is analytic for \z\ < z^{a). The crux of the problem is that 
appears in the summands of (7), and this means that i?(l,0), as it is defined there, is not 
ahsolutely-uniformly convergent for \z\ < z^{a). 

This property can be satisfied for small z (for example, taking \z\ < z^{ca) for a sufficiently 
large constant c), allowing for a resummation^, but unfortunately it is unclear how to extend 
this result to \z\ < z^{a). 

4 Order of the phase transition 

Even though Conjectures 1 and 2 remain unproven at this point, the results of the preceding 
sections are satisfying: we have defined a new model of polymer adsorption, derived functional 
equations satisfied by corresponding generating functions, solved those functional equations, 
and determined the free energy of the model in the limit of infinitely long polymers. This was 
the initial goal of this paper: to find a model more general than directed or partially directed 
walks, which is nevertheless still solvable. This models might then, in some sense, be ‘closer’ to 
SAWs than previously-solved models. We now have the scheme 

fully directed walks C partially directed walks C prudent walks C SAWs. (12) 

We note here that all the existing models discussed in Section 1.1 have second-order ad¬ 
sorption transitions [12,13,25]. Moreover, the general SAW model in two dimensions is also 
conjectured [14] to exhibit a second-order transition. Where it is known exactly, the order of 
these transitions is of second-order for both loops and tails, irrespective of the location of the 
endpoint of the walks. 

^We refer interested readers to [1, pp. 2274-5], where a similar problem was encountered in the enumeration of 
3-sided prudent polygons by area. In that case, however, the generating function was a sum of rational functions 
of a single variable, and is thus much simpler than the problem in this case. 
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Given (12), the nature of the adsorption phase transition for prudent walks (as described in 
Conjectures 1 and 2) may at first seem puzzling. Prudent walks, which undergo a first-order 
transition, are both a subset and a superset of models which undergo second-order transitions. 

A qualititative way of understanding this is to consider the relative location of the singu¬ 
larities for partially directed walks and the two prudent models. In the desorbed state, the 
dominant singularities are given by z\ < z\ < zf^^ = y/2 — 1. On the other hand, in the 
adsorbed state, the equation determining the free energy is given by the solution of the same 
equation for all three models, namely (11). Combining this with the fact that partially directed 
walks undergo a second order transition necessarily implies that the two prudent walk models 
must undergo a first order transition at their respective critical adsorption values. 

The following (rough) conjecture is another approach to understanding this feature. 

Conjecture 3. For a given (infinite) subclass C of self-avoiding walks in an upper half-plane 
with starting point on the surface, let lo G C and define h{u}) to be the maximum height above 
the surface reached by u). Then let (fin) be the average of h{uj) over all walks oj G C of length 
n. If walks in C, equipped with a surface interaction associated with steps along the surface, 
exhibit an adsorption phase transition, then the transition will be first-order if (hn) = Q{n), 
and second-order if {hn) = 0{n^^^). 

The exponent 3/4 is widely expected [19] to be the value v characterising the size of self¬ 
avoiding walks in two dimensions. At present there are no known models with an exponent iz 
strictly between 3/4 and 1, and so we do not attempt to cover this range in the conjecture. 

The intuitive reasoning behind this conjecture is as follows. A first-order phase transition 
characterises a sudden, violent change - a system on one side of the transition point is markedly 
different from one the other. Second-order transitions, while still characterising a change, are 
smoother, and less sudden. 

If a given walk model has {hfi = 0(n), then when desorbed, the bulk of an average walk 
must be far removed from the surface. But when adsorption takes place, a significant portion 
of the walk must now be very close to the surface. This dichotomy between the desorbed and 
adsorbed regimes suggests a sudden, noticeable change must take place at the critical point, 
implying a first-order transition. 

On the other hand, when {hfi = 0{n^/^), the distance between the bulk of a walk and the 
surface is small, relative to the length of the walk. Thus the difference between the adsorbed and 
desorbed phases (at least, close to the critical point) is less dramatic, suggesting a second-order 
transition. 

4.1 Directed models 

In this section we briefly consider Conjecture 3 in the context of existing, directed models. 
We will discuss four such models, all restricted to the upper half-plane, with an impenetrable 
surface in the horizontal line y = 0. 

• NE-directed walks, i.e. walks on the square lattice consisting only of north and east steps; 

• Dyck paths, i.e. walks on a 45° rotated square lattice (or, just as easily, the triangular 
lattice) consisting of north-east and south-east steps; 

• Motzkin paths, i.e. walks on the triangular lattice consisting of north-east, east, and 
south-east steps; and 

• partially directed walks, i.e. walks on the square lattice consisting of north, east and south 
steps. 
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In the literature, the terms Dyck path and Motzkin path are usually used to refer to walks which 
start and end on the surface. To avoid confusion we will instead use the terms loops and tails 
to refer to walks which end on the surface or at any height, respectively. 

In the following table we summarise the results for the asymptotics of the average height 
{hn) of walks of length n, as well as the order of the adsorption phase transition. (Note that 
for models which allow steps parallel to the surface, one can use edge- or vertex-weights on 
the surface, but this does not affect the order of the phase transition.) The calculations for 
NE-directed walks are trivial, and the average heights given there are exact for all n, not just 
in the large n limit. The numerical estimates are ours, and are derived from analysis of long 
(4000 or 5000 terms) series. 


Model 

Loops 

Tails 

NE-directed 

0 

no phase transition 

re/2 

first-order 

Dyck [6,12, 17] 

second-order 

(\/Erlog2) re^/^ “ i + 
second-order 

Motzkin [13,26] 

^1^1/2 _ 3 + 0 (n-V2) 

second-order 

1.418632 re^/2 _ 1,5000 0(re"^/2) 

second-order 

partially directed 
[25] 

1.376996n^/2 _2,91422-KO(n-i/2^ 
second-order 

1.908922 re^/2_2.9i422-HO(re"^/2) 
second-order 


We note that for Motzkin paths and partially directed walks, the ratio of the leading ampli¬ 
tudes of tails and loops seems to be the same as that of Dyck paths, namely 2 log 2 ~ 1.38229. 
This leads us to conjecture that the average height (hn) of Motzkin tails of length n is asymp¬ 
totically 

(hn) = ^2Y^log2^ ^ -h 0(re"^/2). 

If true, this can very likely be obtained in the same manner as the exact results given above. 

More importantly to us, we see that all of these models (except the trivial case of NE- 
directed “loops”, which display no transition at all) obey Conjecture 3. The single model with 
average height 0(n), NE-directed tails, undergoes a first-order adsorption transition, while all 
the others, having average height 0(n^/^), display second-order transitions. 

4.2 Prudent models 

We now return to prudent walks, and investigate the average height of walks of length n and 
how this relates to the order of the adsorption phase transition. We begin with a brief outline 
of how the asymptotics of the average height (hn) can be calculated. 

Let Cn{X) be the number of walks in a class C which reach a maximum height A above the 
surface, and define the bivariate generating function 

C(z;v) = '^Cn{\)z^V^. 

n,A 

Say C{z\ 1) has radius of convergence p. Then 

. EAACn(A) l[z-]C{z-v)U=, 

^ EACn(A) MC(z;l) [z-]C{z-l) ’ 
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where the last equality follows from the uniform convergence of C{z;v) on (at least) {z,v) G 
[Oj/o) X [0,1]. For all the models discussed in this article, we expect to find 

r -r.. d 


[z^]C{z;l) ~ An'^p- 
for constants A, B,'y,y, implying 


and 


|y-l 


BnA p 


V = 1 


{K)=Q{rB 


As per Conjectures 1 and 2, we expect that both tails and loops display first-order transitions. 
For tails, it is easier to consider the average height of the endpoints of walks, (e^), rather than 
the maximum height. Note that since the maximum height of any walk is at least as great its 
endpoint, we have (/i„) > (e„). 

The variable v in W{z\ l,u; 1) tracks the endpoint height of walks. It thus appears that we 
must investigate the asymptotics of 


U = 1 

However, we are fortunate in that W satisfies a set of properties which allow us to sidestep 
this analysis. Theorem IX.9 of [9] states that if a generating function F{z',u) satisfies a certain 
meromorphic schema, then the random variable Xn with probability generating function 


Pn{u) 


[z^]F{z,u) 
[z^]F{z, 1) 


has mean and standard deviation which are asymptotically 0(n). The random variable of 
interest here is Cn, and it has probability generating function 


Pn{v) 


[z^]W{z-,l,v,l) 

[z^]W{z-iAA)' 


We thus need to show that W{z-, l,v; 1) satisfies the three conditions set out in Theorem IX.9 
of [9]. For completeness we state the details of that theorem in the Appendix. 

Condition (i) amounts to showing that W{z] l,u; 1) has an isolated simple pole at {z,v) = 

(^ 1 , 1 )- 

One can factor out the singular term C{z;vA,v)~^ from the generating function and show 
that what remains is analytic on, for example. 


{z,v) G P* = {\z\ < 0.405} X {|u - 1| < 0.01}. 


(Recall that the dominant singularity at u = 1 is a simple pole at z\ ~ 0.403.) The most 
complicated part of this calculation involves showing that the infinite sum (7) is absolutely- 
uniformly convergent on P* for u = 1 and for u = vA, but even this only requires one to derive 
simple bounds for the individual factors of Ti and I. 

Condition (ii) is a straightforward derivative, and Condition (in) can be checked via implicit 
differentiation of C{z;vA,v). It follows that (e^,) = 0(n), and thus the same is true of {hn)- 
For prudent loops, the height is tracked by the variable u in R{z;u,d;l). The dominant 
singularity at rt = 1 is again a simple pole, and so we are still able to use Theorem IX.9 of [9]. 
The singular term (1 — (1 — 2 — -|-A)A(ttA))“^ can be factored out, and what remains can be 

shown to be analytic on, for example, 

{z, m) G = {\z\ < 0.413} X {|ri - 1| < 0.001}. 

(The dominant pole is at zf ~ 0.412.) Conditions (ii) and (in) are again straightforward. We 
thus have (hn) = 0(n). 

We summarise these results in the following lemma, which is given without proof. 
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Lemma 4. The average endpoint height {cn) and total height (hn) of n-step 2-sided prudent 
tails in the upper half-plane are both asymptotically linear in n. The average total height (hn) 
of n-step 2-sided prudent loops in the upper half-plane is also asymptotically linear in n. 

One way of understanding the above lemma is to note that prudent walks, although seem¬ 
ingly undirected, feel the effect of a “pseudo-force”. The step set of prudent walks is unrestricted 
in the sense that the walk is able to step in all four lattice directions, while maintaining the 
prudent condition. However, a 2-sided prudent walk is unable to take the pairs of steps south¬ 
west or west-south. This has the effect of biasing the step set in a direction at an angle to the 
surface. This bias can be interpreted as a “force” acting on the end of the polymer, pulling it 
in the direction of the pseudo-force. 

The effect of a pulling force on a polymer has been considered in the case of self-avoiding 
walks [10] and partially directed walks [20,22]. When the component of the force perpendicular 
to the surface is sufficiently large, desorption by the force occurs and the transition is first-order 
for both models. Furthermore, in the case of partially directed walks, it is known that there 
exists a critical value of the angle, below which adsorption is enhanced, and desorption cannot 
be induced by the force. In this regime, thermal desorption can still occur, and that transition 
is second-order. 

This “pseudo-force” may be one mechanism which leads to the observed first-order transi¬ 
tions for 2-sided prudent walks. 

5 Conclusion 

In this paper we have studied two related models of a polymer interacting with an impenetrable 
surface. The lattice walks considered, namely 2-sided prudent tails and loops, are able to step 
in all four directions on the lattice, unlike many previously-studied directed models like Dyck 
paths and partially directed walks. In order to analyze these walks we have used the generating 
function of 2-sided prudent walks in the upper half-plane. The solution to this generating 
function is obtained by applying the iterated kernel method to a certain functional equation. 
By studying the singularity structure of this generating function, we obtain thermodynamic 
quantities for prudent tails and loops. 

In both cases we locate the critical value of the interaction parameter, and observe a first- 
order adsorption transition. We relate the order of the phase transition in this model with those 
of other models like adsorbing Dyck paths and SAWs, and argue that the difference arises from 
a “pseudo-force” imposed by the restricted step set of prudent walks. 


Appendix 


Theorem (Theorem IX.9 of [9]). Let F{z,u) he a function that is bivariate analytic at (z,u) = 
(0,0) and has non-negative coefficients. Assume that F{z, 1) is meromorphic in z < r with only 
a simple pole at z = p for some positive p < r. Assume also the following conditions. 

(i) Meromorphic perturbation: there exists e > 0 and r > p such that in the domain T> = 
{|z| < r} X {|tt — 1| < e}, the function F{z,u) admits the representation 


F{z,u) 


B{z, u) 
C{z, u) ’ 


where B{z,u), C{z,u) are analytic for {z,u) G V with B{p, 1) ^ 0. (Thus p is a simple 
zero of C{z, 1).) 
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(n) Non-degeneracy: one has dzC{p,l) ■ duC{p,l) / 0, ensuring the existenee of a non¬ 
constant p{u) analytic at u = l, such that C{p{u),u) = 0 and p{l) = p. 


(Hi) Variability: one has 


where ^{g{u)) = g''{l) + g'{l) 



/ 0 , 


Then the random variable Xn with probability generating function 


Pn{u) 


[z^]F{z,u) 

[z^]F{z,l) 


after standardization, converges in distribution to a Gaussian variable, with a speed of conver¬ 
gence that is The mean and standard deviation of Xn are asymptotically linear in 

n. 
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